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Abstract 


We present an autonomous phase-plane describing the evolution of field equations containing 
an agegraphic dark energy in 5D Brans- Dicke cosmology. To observationally verify the numerical 
results, we simultaneously solve the equations by constraining the model parameters with SNe la 
data. We find conditions for the existence and stability of the critical points (states of the universe) 
and numerically examine the cosmological parameters. We also investigate the model by means of 
statefinder diagnostic. 
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1. INTRODUCTION 


Various cosmological observations such as Cosmic Microwave Background (CMB) 


Supernova type la (SNIa) jjl, 0, Weak Leasing 




Baryon Acoustic Oscillations (BAO 


fl. 


2 \, 


2dF Galaxy Redshift Survey (2dFGRS) [7| at low redshift and DEEP2 redshift survey [8| at 
high redshift, have provided cross-checked data to determine cosmological parameters with 
high precision. These parameters imply that our approximately 13.7 year-old universe is 
nearly spatially flat, homogeneous and isotropic at large scale, i.e. a Friedmann-Robertson- 
Walker (FRW) universe with zero curvature, and has entered an accelerating phase since 
z ~ 0.46 { 4 ]. Moreover, according to AC DM model, the universe consists of 0.046 baryonic 
matter, 0.228 non-relativistic unknown matter, namely dark matter (DM), and a remarkable 
amount of 0.726 smoothly distributed dominant dark energy (DE) 2j. The equation of state 
(EoS) of DE, is the main parameter which determines the gravitational effect of DE on the 
evolution of the universe, and can be measured from observations without need to have a 

;hat the EoS of DE lies in a narrow range 

around w 


dehnite model of DE. Strong evidences imply t 

n n 

— 1 and has a smooth evolution [9[, |10| . 


Theoretically, one can classify the 


EoS of DE with respect to the barrier w = —1, namely the phantom divide line (PDF )[ll|. 
That is, DE with the EoS of tc = —1 employs the cosmological constant. A, with a constant 
energy density. The case with dynamical EoS where w > —1, is referred to as quintessence 


12|. and w < —1 corresponds to phantom energy [^, jf^ . 


On the other hand, the problem of DE, its energy density and EoS parameter is still 
an unsolved problem in classical gravity and may be in the context of quantum gravity we 


achieve a more inclusive insight to its properties 


15|. The holographic dark energy (HDE) 


model is an attempt to apply the nature of DE within the framework of quantum gravity 
[3, ill I- The holographic principle states that the number of degrees of freedom describing 
the physics inside a volume (including gravity) is bounded by the area of the boundary which 
encloses this volume and thus related to entropy scales with the enclosing area of the system 
[ 1 ^ . Since the entropy scales like the area rather than the volume, the fundamental degrees 
of freedom describing the system are characterized by an effective quantum held theory in 
a box of size L with one fewer space dimensions and with planck-scale UV cut-off A 181. 

Among all the cosmological models, the scalar-tensor theories have been widely used to 
explain the late time acceleration of the universe and its relation to the HDE 
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One example is Brans-Dicke (BD) theory where the gravitational constant becomes 
time dependent varying as inverse of a time dependent scalar held which couples to grav¬ 
ity with a coupling parameter u. Many of the cosmological problems can be successfully 
explained by using this theory. Alternative 
cosmic acceleration and phantom crossing 


the higher dimensional theories may explain 
4^ . j^. The existence of extra dimensions is 


required in various theories beyond the standard model of particle physics, especially in the¬ 
ories unifying gravity with the other fundamental forces, such as superstring and M theories 


48|-|m|. 


ains the observational data and has been studied 
Kiore recently, a new dark energy model. 


While the successful HDE model exp 
widely by the authors ( see for example 
dubbed agegraphic dark energy (ADE) model, has been proposed by Cai The ADE 
is also related to the holographic principle of quantum gravity and takes into account the 
uncertainty relation of quantum mechanics together with the gravitational effect in general 
relativity. 

Here, in an attempt to integrate both the scalar tensor and higher dimensional theories, 
the 5-dim Brans-Dicke (BD) cosmology is studied with the aim to understand the agegraphic 
nature of dark energy in the model. We perform stability analysis and investigate the 
attractor solutions of the model by utilizing the 2-dimensional phase space of the theory. 
Also, we simultaneously best £t the model parameters with the observational data using 
method. This enables us to End the best htted model parameters for the analysis of 

he model with the experiment. In addition, we 


the critical points and also verihcation of 

use the well-known statehnder parameters 5^ the differentiate among cosmo 
These parameters are used to explore a series of DE cosmological models [57|-[60| 


ogical models. 


2. THE MODEL 

The 5-D Brans-Dicke action is given by, 

S = I + Lm^ , (1) 

where cu is a dimensionless coupling constant which determines the coupling between gravity 
and BD scaler held, R is 5D Ricci scaler, (f){x^) is the BD scalar held and Lm is the 
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Lagrangian of matter field. We assume the metric in 5-D space-time in the form of, 

dr'^ 


ds^ = —dr -I- y) 


+ rdVt^ -I- y)dy^, 


( 2 ) 


1 — kr^ 

where the 4-D space time is assumed to be homogeneous and isotropic (FRW universe). In 
the metric, k is the curvature parameter with k = —1, 0,1 corresponding to open, flat, and 
closed universes, respectively. The scaler held cp and the scale factors a and b, in general 
are functions of t and y. For simplicity and plausibility, we assume the extra dimension is 
cyclic, i. e. the hyp ersurface-orthogonal space-like is a killing vector held in the underlying 
5D space-time Hence, all helds are function of the cosmic time only. Note that, the 
functionality of scale factor b on y, either can be eliminated by transforming to a new extra 
coordinate if 6 is a separable function, and or makes no change in the following equations if 
b is the only held that depends on y. Beside, in the compacting extra dimension scenarios, 
all held are Fourier-expanding around y, and henceforth one can have terms independent of 
y to be observable, i. e. physics would thus be ehectively independent of compactihed hfth 


dimension 




Variation of the action ([T]) with respect to metric ([2]) yields the following held equation. 


4a; 


k 


bH 


1 b 


P 


- — (j) I H -r^H H- - I — —(j) -\- — — H(j)(j) -|- —— -r^Hcpcj) — — PBj 


2oj 


+ + t + V(2i/ + 

4a; V a bH ba^. 


2uj 


2u bH 
b 

m' 


2uj 


2\2 


= Pb, 


-^ch^(2H^ + H+X)-^^ 


4a; 


4> + ‘iHcp -\- 


Au 




bH 




2a; V a 


4 V u 

+ ^H^ + 


O tt2 ^ nTT2 b 




3b 


where H = -. In the orthonormal basis = dt,e^ = adx'^ and = bdy, the stress-energy 


= 0 , 


5 _ 


(3) 

(4) 

(5) 

( 6 ) 


tensor can be considered as. 


T ’A rpA I I rpA I 

B — ^ B \ hulk “r _B \ brane 

where T^\i,uik is the energy momentum tensor of the bulk matter and 

TB\buik = diag{-pB,PB,PBHB,qB) 

The second term corresponds to the matter content in the brane (y=0) 

Tilbrane = diag{-p,p,p,p, 0) 


(7) 

( 8 ) 

(9) 
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Assume that the 05 component of the energy-momentum tensor vanishes, which means that 
there is no flow of matter along the fifth dimension. Therefore the nonzero elements of the 
5D stress-energy tensor are 


rp , ^i.y) 

Tqq — Pb + ^ p. 

(10) 

rp 2 1 2 ^iy) 

Tii — a Pb + a , P- 
b 

(11) 

T 55 = b^QB- 

(12) 

We assume that the BD scalar field 0, is in power law of the scale factor a{t), in the form 

of 0 cx a". In the next section, we apply the ADE model in 5-D Brans-Dicke theory. The 

conservation equations for the dark energy and matter field in the universe are respectively. 

Pb + ^H{pb + Pb) + — 0, 

(13) 

pm “1“ pjyi 0 , 

(14) 


The agegraphic dark energy model with the dark energy density is given by 

Pa = 




2^2 ’ 

where T is chosen to be the age of our universe and given by 

T= f dt = 


(15) 


da 

lU' 


(16) 


In the framework of Brans-Dicke cosmology, we write the agegraphic energy density of the 
quantum fluctuations in the universe as 


PA = 




(17) 


3. STABILITY ANALYSIS-PERTURBATION AND PHASE SPACE 


The structure of the dynamical system can be studied via phase plane analysis, by intro¬ 
ducing the following dimensionless variables. 




, hlyV TT')rr')) Jjn n 1 1 T T ’’ 


Aujpm. „ n 






hH' 


(18) 
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The field equations 


in terms of the new dynamical variables become, 


~ ~^m(3 + 2a + rife + 2-^) 

Sl'^ = + -^) 

n 

~ ~2r2fc(l + -^) 

n; = ^ - 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 


where prime ” / ” means derivative with respect to ln{a). In the above equations we have 


H 

IP 


1 + - 


1 ^ a^(2 + n;) 


+ a + 


By using the Friedman constraint in terms of the new dynamical variables: 

fyfjJ 2 a ' 

— VLm — + 2a(l-+ 126(1 + —) = 1 


(23) 


(24) 


the Eq. ffT^ - fl22l) reduce to, 


Vl’ = -Vtr 


3 + 2a + fib - 


1 + 


1 + 


a^(2 + u) 


acu. 


+ a + 12m + 12 a ~ 2a(l- —) 


— 126(1 + 
12'a = -212a 


2a\ 

y) 

\/iXv 

n 

- S2t(l+^) 


1 / a^{2 + oj) ^ ^ ^ . auj, 

1 H-^-h a + 12m + 12 a — 2a(l- 


1 + - 


(25) 


(26) 


= 3 (- 1 -^ + ^^-126)12a + 2 + ^ 


1 + - 
^3 


3 n 
., 2 /' 


1 + - 
^3 


1 a^{2 + u) ^ ^ ^ _ aw, ^ _ 2a, 

1 H- - -h a + 12m + 12 a — 2a(l- —) — 126(1 H —— 


2(3^ / 

+ 2126(1 H—^) + 2(2 + 126)a + (4 + oj)cP — 126^126 


1 + 


1 + 


a^(2 + w) 


aw. 


2a, 


+ a + 12m + 12 a — 2a(l- —) — 126(1 + ^ 


(27) 


It is more convenient to investigate the properties of the dynamical equations (l2^ - fl27)) than 
equations Eq. flT^ - fl2^ . In stability analysis, the above equations can be solved to find fixed 
points (critical points). These points are always exact constant solutions in the context of 
autonomous dynamical systems. The critical points in this model are highly nonlinear and 
depend on the stability parameters. In addition, in stability analysis the expressions for the 
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critical points and eigenvalues are long and cumbersome such that the usual procedure for 
stability analysis is not possible. To overcome the problem, we solve the above equations by 
best htting the stability and model parameters and initial conditions with the observational 
data for distance modulus using the method. This helps us to find the solutions for the 
above equations and conditions for the stability of the critical points that are physically more 
meaningful and observationally more favored. With the simultaneous stability analysis and 
best fitting, we hnd two critical points. The best htted stability and model parameters and 
initial conditions are a = 0.6, u = —2.9,n = 1.2, h2A(0) = .7, h2fe(0) = .6 and f2m(0) = .27. 
The properties of the two best fitted critical points are given in Table I. 


TABLE I: Best fitted critical points 


CP 

(^A? 

stability 

CPI 

(0,0, 3.5) 

Stable 

CP2 

(0,4.7,0.9) 

Unstable 

CP3 

(0,0,-0.6) 

SaddlePoint 

CP4 

(0,16.3,6.7) 

SaddlePoint 

CP5 

(0.2, 0.2,-0.5) 

SaddlePoint 

CP6 

(0.9 + 27,0.2,-0.5 + 0.91) 

SaddlePoint 

CP7 

(0.9 -27,0.2,-0.5 - 0.97) 

SaddlePoint 

CPS 

(3.6,-2.8,-1) 

Unstable 

CP9 

(4.8,-3,-0.5) 

Unstable 


In Fig. 1, the trajectories leaving the unstable critical point CPI in the past in the phase 
plane is shown going towards the stable critical point CP2 in the future. The best htted 
model parameter trajectories is also shown by green dashed trajectory. 
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Fig. l:The 3-dim and 2-dim phase plane corresponding to the critical point 


4. COSMOLOGICAL TESTS AND STATEFINDERS 


In order to understand the behavior of the universe and its dynamics we need to study the 
cosmological parameters. We have best htted our model with the current observational data 
by the distance modulus test. The cosmological parameters analytically and/or numerically 
have been investigated by many authors for variety of cosmological models. Simultaneously 
best htting the model with the observational data gives us a better understanding of the 
solutions and the dynamics of these parameters. 

We begin with the agegraphic energy density given by equation (l28ll . Taking derivative of 
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(|28|) with respect to cosmic time and substituting into conservation equation (fT3D . in terms 
of new dynamical variables we obtain 

2 a 2^/ni 

1 o o ^b‘ 

6 6na 

With the best htted stability and model parameters, the EoS parameter for agegraphic 
dark energy can be obtained. One has to know that the contribution of extra dimension 
into the formalism is via the presence of parameter a in the held equation for Oa- 


I I I Mllllllllll lllil^M 
/ j I 

III I I I I IIIIIIIIIIIIIIIIII^^^^^^^H 



l.xlO-‘ 1.x 10“ 1.x 10' 1.x 10^ 1.x 10“ 1.x lO** 

Z 


Fig. 2: The best htted EoS parameter for agegraphic dark energy, ica, plotted as functions of redshift 
Among cosmological parameters, the ehective EoS parameter and statehnders are given 


by 


w. 


eff 


= -1 


3(1 + 1) 


(1 + +(+:+ + a + - 2a(l - - ni.(l + y)] (29) 


— are discussed here, where g in r and s is the deceleration 


, r = - 3g - 2 and s = 

parameter and ^ in r in terms of new dynamical variables for exponential and power law 
cases can be obtained by taking derivative of H. 


In Table II, the best htted values of the ehective EoS and statehnder parameters of the 
two critical points are given. 
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TABLE II: Best fitted values of effective EoS and statefinder parameters 


parameter 

Weff 

q 

r 

s 

LOA 

Current Value 

-.85 

-.8 

-3.2 

1.065 

-1.52 

CPI 

-.74 

-.61 

-2.4 

1.065 

-.97 

CP2 

0.46 

1.75 

3.7 

1.415 

.67 


Fig. 3 shows the dynamics of the best htted effective EoS parameter against redshift z. 
From graph or Table II, the best htted current value of EoS tCe// is —.85. The graph also 
shows that universe undergoes phantom crossing twice in future. From a stability point 
of view it shows that the universe starts from an unstable state in the higher redshift in 
the far past and tends to a stable big rip state when the the scale factor of the universe 
becomes inhnite at a hnite time in the future. 


1 

0.8 

0.6 

0.4 

0.2 

0 

-0.2 

-0.4 

- 0.6 

- 0.8 


Fig. 3: The best htted ehective EoS parameter tCe//, plotted as 
functions of redshift,in the solid line) presence and dotted hne)absence of the extra dimension. 

Fig 4 shows the best-htted trajectories of the statehnder diagrams {s,g} and {r, s}. 

From {s, q} graph it can be seen that the best-htted trajectory is currently between SCDM 
and SS state with {s,q} = 1.065,-0.8. Also from {r, s} graph we see that the universe 
passed the LCDM. 



10 -' 10 “ 10 ‘ 10 ^ 10 ^ 10 '* 
-■z 
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Fig. 4:The best fitted statefinders parameters {s — q} and {r — s} 


In Fig 5 we depict the corresponding dynamical behavior of the statefinder {r, s} against 
N = ln{a). From Figs. 4 and 5 we observe that the universe starts its journey from 
unstable state in the past, passed the current state and eventually reaches a stable state in 
the future. Note that in Fig 4, the extreme points of the statefinder trajectories corresponds 
to the state that the universe crosses phantom divide line in future. 



1.765-, 

1.415-' 

1.065 

0.715- 

0.365- 

0.015- 

-0.335- 

-0.685- 


unstable point 


. \stable point \ 


current value 


-5 -4 -3 -2 -1 0 1 2 3 4 5 

N 


Fig. 5:The best fitted statefinder parameters r and s plotted as function of = ln{a). 
fig 6 shows the dynamical behavior of the and for best htted both parameters 

and initial conditions 
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Fig. 6: dynamical behavior of and hlfc plotted as functions of redshift 

Fig.8 shows the dynamic of the (we//) against redshift z in the presence and absence 
of the extra dimension. As can be seen from Fig. 7, there is a slight difference between the 
trajectories of these two cases, we ca see from £g 7 that the current value of EoS parameter 
in absence of extra dimension is {(^eff) — —0.65 while in presence of extra dimension 
~ —0.85 which is more consistent with observational data 



Fig.8: dynamical behavior of effective EoS parameter {uieff) plotted as functions of redshift 
in the solid line); presence of extra dimension and dash line): absence of extra dimension 

Fig.9 shows The distance modulus ^{z) plotted as function of redshift for the best fitted 
parameters 
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The distance modulus fi{z) plotted as function of redshift for the best htted parameters 


5. SUMMARY AND DISCUSSION 

We have constructed an autonomous phase-plane in 5D Brans-Dicke cosmology in the 
presence of agegraphic dark energy and cold dark matter. We have improved the phase-plane 
analysis by constraining the stability and model parameters by SNe la observational data for 
distance modulus using method. In a dark energy dominated universe and with the phase 
space and stability analysis, we hnd two critical points, and the best htted trajectory from 
unstable critical point to the attractor. A quantitative analysis of the dynamical variables 
and physical parameters in the theory is presented. While both effective and ADE EoS 
parameters depend on the BD scalar held, model parameter n and the extra dimension 
b, the best htted trajectories are shown in Fig. 7 in the presence and absence of extra 
dimension, we see that in presence of extra dimension the dynamical evolution of Wg// is 
more consistent with observational data. The results show that the universe has passed 
through SCDM and LSDM phase and approaching the SS phase in future. 
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